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Abstract 



How light can the stop be given current experimental constraints? Can it still be 
lighter than the top? In this paper, we study this and related questions in the context 
of gauge-mediated supersymmetry breaking, where a stop NLSP decays into a W, b and 
gravitino. Focusing on the case of prompt decays, we simulate several existing Tevatron 
and LHC analyses that would be sensitive to this scenario, and find that they allow the 
stop to be as light as 150 GeV, mostly due to the large top production background. 
With more data, the existing LHC analyses will be able to push the limit up to at least 
180 GeV. We hope this work will motivate more dedicated experimental searches for this 
simple scenario, in which, for most purposes, the only free parameters are the stop mass 
and lifetime. 
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1 Introduction 

In supersymmetric extensions of the Standard Model, there are typically many reasons to 
suppose that the stop is the lightest squark. These reasons include: electroweak-scale baryo- 
genesis in the MSSM (for a recent analysis, see [1]); the little hierarchy problem (see, e.g., [21 [3] 
for a recent discussion of this); and the fact that light stops arise naturally in the MSSM - 
renormalization-group running and level-splitting generally pushes the third generation lighter 
than the first and second generations (see |1] for more details). Finally, probably the most 
important motivation of all in the LHC-discovery era - light third generations can lead to very 
interesting and less explored collider signatures. 
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Figure 1: The NLO+NLL stop pair production cross section at the Tevatron (left) and 7 TeV 
LHC (right) as a function of the stop mass. The values of ti cross sections are indicated as 
well. For more details, see appendix IB.ll 

Light stops in theories of gauge-mediated supersymmetry breaking (GMSB) are an espe- 
cially interesting and motivated possibility. As is well known, gauge mediation is an appealing 
supersymmetric scenario: it automatically solves the flavor problem, and it generates phe- 
nomenologically viable soft masses. In such theories, the lightest superpartner (LSP) is always 
a nearly-massless gravitino G. Assuming i?-parity, the next-to-lightest superpartner (NLSP) 
decays in a universal fashion to the gravitino plus its Standard Model partner. Recently, a 
model-independent framework for general gauge mediation (GGM) was established in |^j. 
In GGM, essentially any MSSM superpartner can be the NLSP. So it is interesting to consider 
the case that the NLSP is the lightest stop i. The dominant decay of the stop in such a 
scenario is 



Intriguingly, despite the fact that this possibility has been known for more than a decade (TJ 
m [9] , no searches have addressed it explicitly. And this scenario is far from being obviously 
excluded. 

In this paper, we will focus on the following simple question: how light can the stop NLSP 
be without being in conflict with existing data? In particular, can the stop be lighter than the 
top? Since the stop is colored, stop-antistop pairs have sizeable production cross sections at 
hadron colliders, especially if the stop is light. Still, they can be missed if their decay products 
have a large Standard Model background. Indeed, ti production (where t — )■ W~^b) has a very 
similar signature to tt* production, with a much larger cross section (see figure [1]). Meanwhile 
the uncertainties on the top cross section, both experimental and theoretical, are of the order 
of 10%. As a result, the stop signal may not stand out in ti cross section measurements that 
use simple cuts and event counting. On the other hand, more sophisticated measurements of 
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Figure 2: Diagrams contributing to the stop NLSP decay in GMSB in the simphfied scenario 
where all the other superpartners are heavy. 

the ti cross section and other properties of the top, especially those that require the detailed 
reconstruction of the ti event, may dismiss the stop events as background. One of the main 
goals of this paper is to re-analyze some of the existing Tevatron and LHC measurements of 
tt-like samples and find out to what extent they constrain the stop NLSP scenario. 

For simplicity, we will study a minimal spectrum, consisting solely of stop NLSP t and a 
nearly massless gravitino G. We will take all the other states of the MSSM to be decoupled, 
which is technically possible in the GGM parameter space. This results in two simplifications. 
First, we neglect the contribution from similar decays of the second stop mass eigenstate or the 
production of stops from the decays of other colored states (such as gluinos). Therefore, for 
scenarios in which these additional particles are relatively light, the stop NLSP limits derived 
here should be thought of as somewhat conservative. Second, the decoupling assumption fixes 
the diagrams contributing to the stop decay process (11.11) to those shown in figure O This 
will not limit the generality of our conclusions since the basic kinematic properties of the stop 
events do not depend much on the relative contributions of the various diagrams |8]. 

We start in section [2] by discussing the properties of the stop NLSP decay process. In 
section [3] we discuss the existing Tevatron and LHC analyses that may be sensitive to this 
scenario and use several of them for deriving our constraints in section |H Section [S] discusses 
several other types of measurements that may be relevant to stop NLSPs. We conclude in 
section [61 



2 Phenomenology of stop NLSP 
2.1 The decay process 

In the limit that all of the other superpartners besides the lightest stop are decoupled, only 
the diagrams shown in figure [2] contribute to the stop decay process (11.11) . These two diagrams 
arise from the same gauge- covariant derivative, as explained in more detail in appendix |X1 
The resulting matrix element is 
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Figure 3: Distributions of the Wh invariant mass for stops with masses (in GeV) 120 (blue), 
150 (black), 172 (pink), 180 (green) and 200 (red). We assume = 173 GeV. 



where yF is the SUSY breaking scale, g is the SU{2) gauge coupling, and Cf, describe 
the mixing in the stop sector. For < rrit, the contributions of the two diagrams are 
comparable, while for higher masses the diagram involving the top starts to dominate and 
eventually reduces to the 2-body decay t — )■ tG. The transition between the 3-body decay and 
the 2-body decay is demonstrated in figure El 

From fl2.ip . we see that our model depends on just two parameters: the mass of the lighter 
stop and the stop mixing angle. For a more general spectrum, diagrams with virtual charginos 
or sbottoms would also contribute, but as was noticed in [8], where a much larger parameter 
space has been explored, the kinematic distributions, such as the invariant masses and 
mbw, do not depend strongly on the assumptions about the spectrum or the stop mixing angle 
(which we will set to = —0.8). We therefore believe that the simplified scenario we consider 
is a good representative of the whole class of stop NLSP scenarios. 

In this paper we will focus on the situation in which the stop decays promptly. More 
generally, the lifetime of the stop is dominantly a function of its mass, the SUSY-breaking 
scale \/F (equivalently, the gravitino mass), and various Standard Model parameters. An 
approximate analytical expression for the stop decay rate for < is [8] 
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sin^ Ow 1287r2m2 F2 



(2.2) 



while for > rrit the decay process gradually starts being dominated by t — (with a 
subsequent t — j- W^b decay) which has the rate 

r = ^ (1 — ^ 1 (2.3) 
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Contours of constant stop lifetime are shown in figure HI We see that, as is generally the case 
in gauge mediation, the lifetime of the NLSP can range from prompt (corresponding to lower 
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Figure 4: Contours of constant cr for stop NLSP decay, as a function of the SUSY breaking 
scale \/F and the stop mass m^. 

SUSY-breaking scales and/or heavier stops) to detector-stable (higher SUSY-breaking scales 
and/or lighter stops). For < m^, prompt decay of the stop requires the SUSY breaking 
scale a/F to be as small as it can possibly be, on the order of 10 TeV. 

It is also important to consider longer lived stops, but we will not do so in detail in 
this paper. Stops that are sufficiently stable that they travel fully through the detector are 
constrained by searches for stable charged or colored particles. The current best limit on 
detector-stable stops comes from ATLAS and corresponds to > 300 GeV [10]. (Very 
long-lived stop NLSPs may also have important consequences for BBN [TTl fT2] .) Even more 
interesting is the intermediate case of a stop that decays at a displaced vertex. As far as 
we know, there are currently no limits on this scenario. This would give rise to signatures 
involving displaced jets and leptons. This could pose interesting challenges for triggering and 
reconstruction, as was recently discussed in a related context in [13]. 

2.2 Kinematic distributions 

In figures [SHHl we plot the distributions of various kinematic quantities characterizing pair- 
production of stop NLSP. These plots were made using a combination of our own code for the 
decay of t — WbG according to the matrix element (12.11) : and Pythia for everything elseQ 
Unless stated explicitly, these distributions are purely parton-level, i.e. they do not include 
showering, hadronization, detector simulation, or any cuts. 

The distributions for the Tevatron and the 7 TeV LHC turn out to be very similar, both 
quahtatively and quantitatively (this happens because the tops and stops are produced rela- 
tively close to threshold), so we present them only for the LHC case. For comparison, we have 

^Note that decaying the W through Pythia has the effect of averaging over the W polarization. This is 
not a problem since the W polarization fractions for the stop are very model dependent anyway, while their 
effect on the basic kinematic distributions is relatively small |5]. 
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Figure 5: Parton-level distributions of the leptons pt (top left), 6-jets pt (top right) and I^t 
in the different channels (bottom) for ti (thick black line) and stops with masses (in GeV) 120 
(blue), 150 (black), 180 (green) and 200 (red). 
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Figure 6: Left: parton-level distributions of m^, eq. (12. 4p . in the lepton+jets channel. Right: 
same distributions after showering (the full ti sample) and applying geometric acceptance, 
lepton selection and the dilepton veto (as defined in [H]). 
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Figure 7: Parton- level distributions of Ht-, eq. ( I2.5p . for the lepton+jets (left) and dilepton 
channel (right). 
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Figure 8: Parton- level distributions of m^^ where the lepton and the h quark come from the 
same top or stop (left), the opposite combination (middle) and random pairing (right). 

also included distributions of the same quantities for ti. 

The transverse mass m\. = —Pr of the W , is determined from the measured quantities 
in the lepton+jets channel as would be appropriate if the neutrino were the only invisible 
particle, namely 

rriT = ^J2{p't^t-Vt-^t) = ^2p^r^r (l - cos A0,,^^) (2.4) 

is defined a^ 

Ht= ^Pt + ^Pt + ^t (2.5) 

leptons jets 

^Note though that the definition of Ht varies between the different experimental analyses, sometimes 
excluding I^t or including only jets. The definition (j2.5p which we use in figure [7] is often also denoted TOoff- 
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Some of the distributions are shown separately for the different decay possibihties of the two 
ly's: dilepton, lepton+jets and all-hadronic (by lepton, we mean e or /i). Note that unhke 
for tt, stop events have significant missing transverse energy (MET) carried by the gravitinos 
even when both W^s decay hadronically, which results in a jets+MET signature. 

Many of these kinematic quantities are used in the experimental analyses that we will 
study in this paper, either to select for ti events, or to distinguish such events from events 
with stops or other top partners. One major exception is mih, the invariant mass of the lepton 
and the 6-quark coming from top or stop decay. The distributions are not used in any 
of the existing analyses, even though they can be extremely useful in distinguishing between 
stop and top events [8] . We will return to met in section 15.11 

One of the most significant differences relative to the top for < rrit is that the b jets are 
much softer (figure [5]). This happens due to purely kinematic reasons: in the decay t — )■ W~^b, 
the b momentum (in the top rest frame) is 

m ( ^ 



P^ = -^(^--^ (2.6) 
2 \ mt 



(where for simplicity we neglected rrih), which creates a Jacobian peak in pf,^T (which gets 
somewhat smeared, as seen in figure [5|, because the tops are produced with finite px)- On 
the other hand, in the case of the 3-body decay t — W~^bG the momenta of the b quarks are 
distributed within 



Pb<^{l-^] (2.7) 



2 



2 V m? 



Additional differences are a significant high-mr tail (in the lepton+jets channel) and the 
somewhat lower values of Ht in the lepton+jets and dilepton channels. In section H] we will 
see how these differences determine the acceptance of the stops in the various existing analyses. 



3 Overview of relevant Tevatron and LHC analyses 

In this section, we will give a brief but comprehensive overview of the Tevatron and LHC 
analyses which are relevant for constraining the stop NLSP scenario. We have divided up 
the discussion into two categories: Standard Model analyses of ti production, where stop pair 
production could be "hiding"; and explicit searches for new physics with stop-like final states, 
for which ti is a major background. 



3.1 tt analyses 

Because of the general similarity between top and stop NLSP signatures, measurements of the 
ti cross section have the potential to also constrain stop pair production. Experimental analy- 
ses of ti production fall into four categories, based on how the two W's decay: dileptonic (5%), 
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Analysis 


Leptons 


6 tags 


Luminosity 


Excluded 

stop masses 




CDF [15] 


2 


> 

> 1 


4.47 fb-i 
4.47 fb-i 


- 
- 


/ 


CDF [16] 


2 


> 

> 1 


5.1 fb^i 
4.8 fb-^ 


- 




DO [17J 


2 


> 


5.4 fb"^ 


not simulated 




CDF [18] 


1 


> 1 


4.3 fb-i 






DO pJ] 


1 


> 1 


5.3 fb-^ 


not simulated 




ATLAS [20] 


2 


> 

> 1 


35 pb-i 
35 pb-i 


- 


/ 


CMS [21] 


2 


Al Al 


36 pb-i 
36 pb-i 


not simulated 
not simulated 




ATLAS [221 [23] 


1 


> 

> 1 


35 pb-i 
35 pb-i 







Table 1: The most recent Tevatron and LHC ti cross section measurements (with cut-and- 
count analyses). The dashes indicate that no ranges of stop NLSP masses were found to be 
excluded at 95% CL by our reproduction of the analyses. The last column indicates the most 
sensitive analyses whose results will be presented in more detail. 

lepton+jets (35%), jets+MET (10%), or all hadronic (50%)! The dilepton channel is gener- 
ally the cleanest: after requiring two leptons, two jets, and MET, the signal-to-background 
ratio for events is typically ~ 2, and after requiring at least one of the jets to be 6-tagged 
it increases to ~ 10. However, due to the small branching ratio of this mode there is sta- 
tistical uncertainty roughly of the same size as the systematic uncertainty (this is true for 
both 5 fb"^ at the Tevatron and 35 pb~^ at the LHC, for both pre-tag and 6-tagged samples). 
The lepton-|-jets channel has a much larger branching ratio, however the pre-tag sample is 
dominated by the W^+jets background, so usually only the 6-tagged sample is relevant (unless 
more sophisticated techniques that assume detailed information about the signal properties 
are used, but those would have less sensitivity to new physics). After a typical set of require- 
ments (a lepton, > 4 jets, MET) the signal-to-background ratio is ~ 4 and the systematic 
uncertainty dominates over the statistical uncertainty (by a factor of ~ 2). The most recent 
Tevatron and LHC ti cross section measurements in the dilepton and lepton+jets final states 
using cut-and-count methods are listed in table [TJ 

■^By lepton we mean e or ^ (including the case of a leptonically decaying t). The jets+MET category 
consists of zero-lepton events where at least one W decays to tv and the t subsequently decays hadronically. 
The all-hadronic category, by contrast, consists of events with no intrinsic MET, where both VF's have decayed 
hadronically. 
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To get a rough idea of the sensitivity of ti cross section measurements to stop pair pro- 
duction in these channels, we can consider the ratio of stop and top production cross sections, 
which can be read off figure [TJ For example, for mi = nit, ^ihl^tt = 0-10 the Tevatron 
and 0.18 at the LHC. Since tt is a dominant background to stop pair production, this gives 
a crude estimate for the signal-to-background ratios, both with and without h tagging. Of 
course, these numbers will be affected, sometimes very significantly, by the differences in the 
kinematic properties between stop and top events discussed in section [2J 

Decay modes without leptons seem less useful because of the large QCD background. A 
possible exception is the jets-|-MET channel. Because of the MET carried by the gravitinos, 
most of the stop events fall in this category. By contrast, the ti background is significantly 
reduced (branching ratio ~ 10% as indicated above). The QCD background (where MET can 
arise due to mismeasured jets) is still an obstacle for a cut-and-count analysis, but a dedicated 
search using for example neural networks may work (similarly to the top mass measurement 
in this channel |2l])- However, designing such a search and examining its feasibility is beyond 
the scope of the present paper. 

Various properties of the top quark have been measured besides the cross section. In par- 
ticular, the top mass has been measured at the Tevatron with the great precision of 1 GeV. 
However, these measurements are not necessarily more sensitive to the presence of stop events. 
Many of them use matrix elements or neural networks for background rejection and are there- 
fore likely to discriminate against stop events. Others, in order to extract the desired quantity 
and/or obtain a pure sample of tops, require a detailed reconstruction of the ti event from the 
observed objects. Since the momentum carried by the gravitinos will be unaccounted for, stop 
events will not reconstruct in a meaningful way and are likely to be rejected as background 
or make a relatively smooth contribution to the measured distributions. These analyses also 
present a practical difficulty for us, since without having access to the code used in the experi- 
mental study, reproducing the behavior of complex algorithms for events which they were not 
designed to treat would be risky. For these reasons, we do not include these measurements in 
our study. 

3.2 Searches for new physics 

Several experimental studies have looked specifically for light stops, but not in the GMSB 
scenario with stop NLSP that we consider here. However, in some cases the stops also had 
tt-like signatures and such analyses can have good acceptance to GMSB stops as well (despite 
the fact that they were optimized for other scenarios). We have listed the relevant searches in 
table El and will now proceed to describe them in more detail. 

At the Tevatron, CDF [25] and DO [261 EZl EH] have considered the situation in which the 
stop decays via a virtual chargino as 

t b£+i> (3.1) 

where the massive sneutrino u can be the LSP or decay invisibly as — )■ i^Xi or ^G. The 
signature here is similar to that of the dilepton channel of (11. ip . The DO search [28] (with 
5.4 fb~^ of data) used distributions of composite discriminant variables optimized for the decay 
(13. ip as a function of the stop-sneutrino mass difference, and excluded stops with masses as 
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I\iL(Xiy bib 


1 OT~\"f" r^Tl c 

IjCjJ LUlio 




I-j UllilliUol 


Excluded 
stop masses 




DO stop search [28] 


2 


> 


5.4 fb"^ 


mi < 123 


/ 


CDF stop search [30] 


2 





2.7 fb"^ 


- 








> 1 


2.7 fb-i 


105 < mi < 120 


/ 




with stop mass reconstruction 


< 150 


/ 


DO stop search [22] 


1 


> 1 


0.9 fb"^ 


not simulated 




ATLAS top partner search [Tl] 


1 


> 


35.3 pb-i 




/ 


ATLAS SUSY search [33] 


> 1 


> 1 


35 pb-i 







Table 2: Recent new physics searches and the stop NLSP masses (in GeV) excluded by them 
at 95% CL according to our analysis. The last column indicates the most constraining analyses 
whose results will be presented in more detail. 

large as 240 GeV for a certain range of u masses (assuming 100% dilepton branching ratio, an 
order of magnitude larger than in our case). 
Another possibility that has been studied is 

i^bxt, xt^W+'^*^Xi (3.2) 

which also gives rise to a ti-like signature. The most recent such study from DO in the 
lepton+jets channel (although with only 0.9 fb~^) did not yield significant exclusion limits [29]. 

Finally, CDF [30| |3T| [32] has considered a scenario which has some overlap with (13. 2p in 
the dileptonic final state: 

i-^bxt, Xt^^^y^i (3.3) 

The CDF study treated the chargino branching ratio to leptons as a free parameter. If the 
chargino decays through a ly as in (13. 2p . then the dilepton channel will have branching ratio 
of 0.11 (as is the case for our stop NLSP scenario). The chargino can also decay through a 
charged higgs, a slepton or a sneutrino, in which case the branching ratio can be as high as 
1. In order to discriminate between stop and top events, CDF has designed an algorithm for 
approximately reconstructing the stop mass. Because the final state contains four undetectable 
particles and the masses of the intermediate chargino and the final neutralino are unknown, 
it is impossible to reconstruct the event rigorously. Nevertheless, after assuming a mass for 
the chargino, the algorithm is able to construct a quantity, "stop mass," whose distribution 
is peaked roughly at the stop mass. For top events, one obtains a much broader distribution 
peaked somewhat above the top mass. Using 2.7 fb~^ of data, CDF were able to exclude 
stop masses of up to 200 GeV for certain chargino and neutralino masses assuming a dilepton 
branching ratio of 1, but only a small range of stop masses around 130 GeV for a branching 
ratio of 0.11. 

Among the searches for new physics released by the LHC experiments to date, two may be 
potentially relevant to light stops. One is the ATLAS SUSY search for events with leptons. 
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6-jets and missing energy [53] • This search assumed the decay channel f l3.2p . It did not have 
sufficient exclusion power for directly produced hght stop pairs (but only for stops produced in 
decays of relatively light gluinos). The search was not optimal for such stops, in part because 
it required the effective mass of the event to be rries > 500 GeV. Another interesting study 
is an ATLAS search for a heavy top partner (a fermion decaying into a top and an invisible 
particle) [H]. This is potentially better suited to the light stop NLSP scenario, because it 
does not impose a hard nics cut. 

4 Constraints from the Tevatron and the LHC 

Using crude detector simulations and object definitions (see appendix O for details), we have 
reproduced a large number of the analyses listed in tables [H and [2] to a reasonable level of 
accuracy. We have then used this to estimate the current experimental constraints on stop 
NLSPs. As can be seen from the tables, most analyses currently set no limit at all on stop 
NLSPs. This is not surprising, given that no analysis has been optimized to this scenario. 
Experimental sensitivity is degraded by a combination of low acceptance, especially for the 
lighter stops, and/or low cross section, especially for the heavier stops. 

In this section, we will describe in more detail those analyses that we have found to be most 
sensitive to the stop NLSP scenario. These are indicated by the check marks in tables [1] and [21 
As we will discuss below, these searches succeed either because they are more accepting to 
soft jets, or because they use more discriminating variables such as mj- or the "reconstructed 
stop mass". 

For the most part, we have focused on the simple "cut-and-count" portions of these anal- 
yses. This is the type of analysis that we can simulate most reliably, and anyway, we expect 
such analyses to be the most receptive to scenarios that differ somewhat from those for which 
they were designed. One exception is the CDF stop search [30], where we were able to repro- 
duce their more sophisticated analysis based on the "stop mass" reconstruction. This actually 
ends up setting the best limit on the stop mass. For the DO stop search [2B] , which also uses a 
more sophisticated approach, we were unfortunately unable to make use of their distributions 
of composite variables without having access to data. Here we will only use the simple cuts 
which precede that procedure (the "selection 1" described in [28]). 

We will discuss the "cut-and-count" analyses in section 14. 1^ and the CDF stop mass re- 
construction in section 14.21 

4.1 Cut-and-count analyses 

For each analysis and each final state, we compute the expected number of stop events using 
the following formula: 



The ffist factor is obtained from our simulation. It is the acceptance for stop pair production 
relative to that for ti. This is shown in figure |9] for the highlighted analyses. The second 
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Analysis 


it* 


bg 


total bg 


data 


CDF ti cross section [15] (pre-tag sample) 


21 


156 


223 ± 16 


215 


DO stop search [28] (selection 1) 


46 


183 


1174 ± 73 


1147 


CDF stop search [30] (fe-tagged sample) 


12.9 


45.0 


52.4 ± 7.2 


57 


ATLAS ti cross section [20] (pre-tag sample) 


17.5 


74.7 


96.3 ±8.6 


105 


ATLAS top partner search [11] 


7.2 


9.2 


17.2 ±2.6 


17 



Table 3: The expected numbers of events for a 150 GeV stop NLSP, tt background, total 
background (and its systematic uncertainty, as discussed in appendix IB. 21) and the observed 
number of events. 

factor is the stop cross section (figure [T]), again relative to the ti cross section. The final 
factor is the expected number of ti events, which we take from the experimental analyses 
themselves, with one important modification - we have normalized the ti numbers of events to 
the NNLOapprox±NNLL cross sections quoted in appendix IB.lt which generally differ slightly 
(typically lower by ~ 10%) from the values assumed in the original experimental analyses. 
Table [3] shows the event yields for the case of a 150 GeV stop. We estimate that the event 
yields are correct to the ~ 10% level (for more details, see appendix [C]). This is sufficient for 
all practical purposes since this is anyway of the size of the theoretical uncertainty on the stop 
cross section. 

Finally, with the expected number of stop events in hand, we estimated the exclusion 
confidence levels (CL) for the stop NLSP as a function of its mass using the frequentist 
method [51]. Our confidence levels include estimates of the systematic uncertainties based on 
the information available in the experimental papers - for more details see appendix IB. 21 The 
95% CL excluded cross sections, relative to the theoretical stop cross section, are summarized 
in figure [TOl For the LHC searches we also present the expected exclusion limits for 300 pb~^ 
and 3 fb~^ of data. We computed them by assuming the number of observed events for those 
luminosities to be the expected number of background events, and the systematic uncertainties 
(in %) to remain the same. This is a conservative, somewhat pessimistic assumption since 
some of the systematics will improve with more data. 

Figures [9] and [10] contain the main results of our paper. Some comments on these results 
are now in order. 

By estimating the number of stop events via the ratios in (14.1 p and the experimental 
prediction for Nfj, many of the systematic errors in our simulations should cancel out. This 
gives us some confidence in the accuracy of our results. Another important check is that the 
raw number of ti events predicted by our simulations generally agrees at the ~ 30% level with 
N^j from the experimental references. (For more details see appendix [Cl) Note that in any 
event, the limits on the stop mass are fairly robust. Changing the acceptance by ~ 10% in 
either direction would not affect the stop mass limit by very much, given the rapid power-law 
dependence of the stop cross section as a function of its mass. 

We find that for very light stops (m^- ~ 120 GeV), the acceptance is affected significantly 
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Figure 9: Acceptances of tt* events relative to those of tt events for the pre-tag sample of the 
CDF tt cross section measurement in the dilepton channel [15] (blue), the DO stop search in 
the e/x channel (up to selection 1) [28] (red), the fe-tagged sample of the CDF stop search in the 
dilepton channel [30] (green), the pre-tag sample of the ATLAS ti cross section measurement 
in the dilepton channel [20] (pink) and the ATLAS top partner search [T3] (gray). 

by the requirements regarding the number and of jets. Recall from figure [5] that the h 
jets coming from light stop decays are very soft. In figure [9], the acceptance is close to 1 only 
for the DO stop search (red curve) because that analysis does not impose a requirement on 
the number of jets. As a result, it gives the best cut-and-count exclusion limit for a 120 GeV 
stop (see figure fTOl) . The acceptance is still relatively high for the 6-tagged sample of the CDF 
stop search (green curve), which requires one jet of 15 GeV and another of 12 GeV. All the 
other analyses in figure [9] (and even more so some of the analyses from tables [1] and [2] which 
are not included in the figure) have stricter jet Ex requirements which result in much lower 
acceptances. This is true for both dilepton and lepton+jets analyses. 

For heavier stops (mj > 150 GeV) the h jets become harder and are able to satisfy the 
selection requirements much more efficiently, but the limits become significantly weaker be- 
cause of the smaller cross sections. The ATLAS top partner search in the lepton-|-jets channel 
has an exceptionally high acceptance relative to tt, because it uses a cut on ttit requiring 
tut > 120 GeV. As shown in figure [6l this eliminates the top background (where rriT peaks 
sharply below rriw since the MET is all from a single W ^ iu decay), while still retaining 
a decent fraction of the signal. This search is still unable to exclude any stop mass range at 
95% CL, but this is due to low statistics. We find that with several hundreds pb~^ of data it 
will have a potential for 95% CL exclusion for stop masses up to about 180 GeV, as shown in 
figure [TOl but beyond that the analysis is limited by (our rather pessimistic extrapolation of) 
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Figure 10: 95% CL excluded cross sections, relative to the theoretical stop cross section. On 
the left we present our limits from the Tevatron analyses: the pre-tag sample of the CDF 
ti cross section measurement in the dilepton channel [I5] (blue), the DO stop search in the 
e/x channel (up to selection 1) [28] (red), and the 6-tagged sample of the CDF stop search in 
the dilepton channel [30] (green). The two thick green lines are obtained from the stop mass 
reconstruction procedure of [30]. The right plot presents our limits from the LHC analyses: 
the pre-tag sample of the ATLAS ti cross section measurement in the dilepton channel [20] 
(pink) and the ATLAS top partner search [Hj (gray), where the solid lines are the actual 
limits (from 35 pb~^ of data) while the dashed and dotted lines are approximate expected 
limits for 300 and 3000 pb~^, respectively. 



the systematic errors. 

4.2 Stop mass reconstruction 

We have implemented the dilepton stop mass reconstruction algorithm of the CDF analysis [30l 
[311 [32] (in a somewhat simplified way described in appendix [D]). Despite the fact that the 
logic of the algorithm uses the assumption that the stop decays as t — > bxt — )■ d^vhy^, we 
find that surprisingly it works also for our stop NLSPs, in the sense of giving a reconstructed 
mass distribution that is much sharper than that of the top and centered at a different value. 
Our reproduction of the distribution from [30] (for mi = 132.5 GeV) is shown in figure [TT](A) 
(for the 6-tagged channel). The same figure shows also the result we obtain for a stop NLSP 
of the same mass. 

We can use the results of the CDF study for setting limits on stop NLSPs by finding pairs 
of NLSP and gravity-mediated stops with similar mass distribution (which would sometimes 
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Figure 11: Normalized distributions of the reconstructed stop mass obtained for the top (sohd 
black), and the gravity- mediated (dashed red) and NLSP (dotted blue) stops from table Hi 
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case 


NLSP 


gravity-mediated 




case 


NLSP 


gravity- mediated 




■mi 








BR 




mi 


mi 




m^o 

Xi 


BR 


A 


132.5 


132.5 


105.8 


47.6 


0.17 




E 


140 


150 


125.8 


58.8 


0.25 


B 


145 


137 


105.8 


47.6 


0.13 




F 


150 


155 


125.8 


58.8 


0.20 


C 


155 


140 


105.8 


47.6 


0.11 




G 


160 


160 


125.8 


58.8 


0.16 


D 


165 


150 


105.8 


47.6 


0.10 




H 


165 


160 


125.8 


58.8 


0.13 



Table 4: Pairs of NLSP-gravity-mediated stops which have equal numbers of events passing 
the selection of the 6-tagged channel of [50] and similar reconstructed mass distributions 
(figure [TTj) . The first table corresponds to figure 2(a) of [30] which assumes chargino mass of 
105.8 GeV and the second table to figure 2(b) with chargino mass of 125.8 GeV. 

require them to have non-equal true masses) and same number of events (which can be obtained 
by tuning the dilepton branching ratio of the gravity-mediated stop). Note that by taking 
this approach we bypass the need to reproduce the analysis of the systematic uncertainties. 
Also, the imperfections in our simulation of the stop NLSP events and the mass reconstruction 
algorithm are likely to be canceled to a large extent by similar imperfections in our simulation 
of the corresponding gravity-mediated stops. 

Table m shows pairs of NLSP and gravity-mediated stops that have same event counts and 
similar reconstructed mass plots (in the fe-tagged channel) which are shown in figure [TTJ With 
these pairs in hand, we can approximately read off the corresponding 95% CL exclusion limits 
from figure 2 of [30] . We include the results in figure [10] (thick green lines, corresponding to 
the two chargino mass hypotheses of [30]). The conclusion is that stop NLSP is excluded for 
mi < 150 GeV. 

Evidently, the more sophisticated approach of the CDF stop search sets the best limit on 
stop NLSPs. This illustrates the power of using more discriminating variables in searching for 
new physics in the ti sample. We are optimistic that with more work the limit could be further 
improved with existing data. For instance, it would be interesting to see what constraints the 
existing DO stop search [2B] can set when going beyond the simple cut-and-count portion 
("selection 1") that we have considered here. 

5 Other types of measurements 

In this section we discuss several additional methods that may be relevant for future searches 
for stop NLSPs. 

5.1 rngb and 6-jet pt 

As has been pointed out already in [S], the stop NLSP has a distinct distribution of the invari- 
ant mass mih which can be useful for reducing the ti background. We show this distribution 
in figure [HI While there exists an ambiguity in pairing each lepton with the b jet that came 
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from the same top or stop, we see that even incorrectly paired cases happen to contribute in a 
similar way to the differences between the distributions. This happens because the lower 
values for the stops can be attributed to a large extent to the lower momenta of the b quarks. 
Unfortunately, it seems that raw data for the rriib distribution of ti samples at the Tevatron 
have not been published since the early measurements with ~ 700 pb~^ |35]. Such data from 
the Tevatron or the LHC may be able to strongly constrain the stop. 

It is also possible that just the analysis of the 6-jet px distribution will have a comparable 
power. As discussed in section 12.21 (see figure E]), the 6-jets coming from stop NLSP decays 
are generally much softer than those from top decays. If some way could be devised to take 
advantage of this separation, this could become a useful way of distinguishing stop events from 
top events. Some potential limiting factors for such an analysis might be: that only jets with 
Pt ^ 12 GeV can be properly reconstructed; that without b tagging ISR jets may contribute 
as well; that b tagging efficiency decreases with pt] and that additional backgrounds become 
important once low-p^ jets are allowed. 

5.2 Displaced decays 

As we mentioned in section |21 the stop NLSP can naturally be long lived. (Another situ- 
ation in which the stop can be long-lived due to a suppressed coupling is the right-handed 
sneutrino LSP scenario [36].) If the stop decays between ~ 100 /im and ~ 0.5 m away from 
the interaction point in the plane transverse to the beam (the precise range depends on the 
particular experiment), the tracks of the charged particles produced from its decay products 
can be identified as emerging from a displaced vertex. For distances above ~ 1 cm, where 
the background of displaced vertices coming from heavy flavor production becomes small, the 
displaced vertices of the stop can even make its signal easier to separate from the ti and other 
standard model backgrounds. Triggering on such events is easy thanks to the presence of 
leptons (from the leptonic decays of the b or the W). Reconstructing the displaced vertices 
is also simple because of the presence of a 6 jet and a lepton or two additional jets, each jet 
typically containing multiple charged particles. 

Even without using the presence of displaced vertices, such events may be included in the 
samples of prompt searches like those we studied in this paper and our analysis will apply. 
But this depends on the details of the logic used in the event selection procedure of each 
experiment. For example, such events may not make it into the sample or may be analyzed 
incorrectly if some part of the procedure uses the assumption that all the relevant tracks in 
the event emerge from the vicinity of the primary vertexjj 

So far, there have not been any dedicated searches for stop decays at displaced vertices. 
However, many of the ideas that would be useful for such a search can be borrowed from 
searches that addressed displaced decays of certain other particles. In particular, the DO 
search for pair produced long-lived neutral particles that each decay into bb within the 
tracker (motivated by the hidden valley scenario) seems relevant. Methods for studying such 
a scenario have been developed also at ATLAS [38], who in addition considered the possibility 

^For example, the DO stop search [2H1 requires the distance of closest approach between the muon track 
and the primary vertex to be < 0.02 cm if the track includes hits in the silicon tracker or < 0.2 cm otherwise. 
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that the decays occur within the calorimeters. The large missing energy present in the stop 
NLSP scenario, as well as the fact that the stops are somewhat heavier than the particles 
considered by DO and ATLAS, may be helpful for further eliminating backgrounds. The DO 
search (with 3.6 fb~^) has excluded signals with cross sections of the order of 10 pb. It is 
therefore plausible that such a method will be efficient for studying the displaced stop NLSP 
scenario. 

Furthermore, unlike in the scenario described in the previous paragraph where the particles 
are neutral, the stops will hadronize with light quarks into either neutral or charged parti- 
cles [7]. In the latter case, they will have tracks characterized by an anomalously large rate of 
energy loss through ionization (dE/dx) which can be measured if they traverse a sufficiently 
long distance within the tracking volume. In fact, this kind of measurement has been already 
used by CMS [39] (with 3.1 pb"^) and ATLAS [10] (with 34 ph'^), although without any 
results relevant to displaced stops. The CMS study did not have the reach to exclude stops of 
any mass. It should also be noted that the study assumed the stops to be stable throughout 
the calorimeters and computed the trigger efficiencies accordingly. If the stop decays earlier, 
the signal should be simulated differently. For example, the missing energy in the scenario 
considered by CMS is carried by the stops themselves (which similarly to muons do not de- 
posit much of their energy in the calorimeters) while in our case it would be carried by the 
gravitinos and sometimes the neutrinos. The ATLAS study required the tracks to match to 
either a reconstructed "muon" in the muon calorimeter or to a cluster in the tile calorimeter. 
As a result, their analysis does not have acceptance for stops that decay before reaching those. 

5.3 Bound states 

An interesting feature of the stop NLSP scenario is that the near-threshold tt* bound state, 
the stoponium, is guaranteed to decay by annihilation since its annihilation rate (into gg) is 

Tannih ~ ^a'^,- ~ (lO'^^ m)"' (5.1) 

oi 

where we evaluated the decay length for m^- ~ mt (compare with figure H]). Such a stoponium 
will be observable at the LHC within a few years as a narrow diphoton resonance from 

{it*) ^ 77 (5.2) 

at invariant mass of slightly below 2m^ (for recent discussions, see [iQlHT]). The observation 
of this signal will allow a precise measurement of the stop mass. The size of the signal and 
the angular distribution will also help confirm the identity and the properties of the stop. 

5.4 Flavor-violating decays and same-sign dilepton signals 

If the MSSM has flavor violation beyond that of the Standard Model (which is not the case 
in gauge mediation), the process 

i^cG (5.3) 
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may become significant and even dominate over f 1 1.11) [7] . While we are unaware of any existing 
searclies of (15.31) . a similar process 

t^cx? (5.4) 

(with a massive neutralino Xi) has been considered by DO |12] and CDF [13]. The CDF 
search was able to exclude stops up to 180 GeV in the case that this process dominates and 

— m^o > 50 GeV. The limit would probably be even stronger in the case of (15. 3p - since the 
gravitino is massless, the c-jet px and the missing energy would both be larger. Since we are 
mostly interested in light stops {rrif < rrit), the contribution of (15.31) can be considered small. 

However, even a relatively small extra flavor violation, for which (11.11) still dominates, 
may lead to the conversion of the stops into antistops when they bind with light quarks into 
neutral mesinos before decaying [7j. This would lead to same-sign dilepton events with jets and 
missing energy which have very little background. The non-observation of such events [HI US] 
will provide support to the gauge mediation scenario if a light stop is observed through (11.11) . 

6 Conclusions 

We hope we have convinced the reader that there is a viable possibility for the stop to be 
light, and hiding in the tt sample. On the other hand, we have shown that some existing and 
ongoing Tevatron and LHC searches do have sensitivity to this scenario. Now is therefore the 
right time to study it and perhaps discover supersymmetry in our backyard while the LHC is 
pushing the other squarks and the gluino to higher and higher masses. 

Besides determining the fate of the light stop NLSP, the relevant searches would also be 
useful for gaining more confidence in the purity of the top sample or maybe discovering a 
different new physics contribution within it. This seems especially important in view of the 
anomalous forward-backward asymmetry observed at the Tevatron. 

Our limit > 150 GeV was obtained by analyzing measurements that were designed 
for other purposes. We have also seen that even without any further optimization, existing 
LHC searches should have sensitivity to the entire range of < rrit with just 300 pb~^. 
These results can almost certainly be improved by performing more dedicated searches for 
stop NLSP at the Tevatron and LHC. For example, our simulation of the simple cuts stage 
of the DO search [28] has already allowed us to set the limit > 123 GeV, but without 
having access to the data we were unable to utilize the full strength of the DO approach 
which used distributions of a set of powerful discriminating variables optimized for a different 
stop scenario. Adapting this kind of search for the stop NLSP scenario may improve our 

^ 150 GeV limit. Even more simply, the signal-to-background ratios for many analyses 
will benefit from allowing for softer jets and/or using more discriminating variables such as 
rriT or m£b to reduce backgrounds. Clearly, dedicated searches for the stop NLSP scenario 
would be very beneficial. 

In this paper, we have focused on prompt decays of stop NLSPs. This has allowed us to 
cast our results in a ID parameter space consisting solely of the stop mass. However, as we 
have discussed, it is natural for the stop to be long-lived. Therefore, a larger 2D parameter 
space is relevant, consisting of the stop mass and lifetime. Most of this parameter space is 
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unexplored territory. For instance, while the case of a detector-stable stop has been analyzed 
in [TU] , we are not aware of any searches which constrain stops with displaced decays inside the 
detector. Furthermore, even in searches optimized for either of the two extreme cases (prompt 
and detector-stable), it is not known how the acceptance degrades as one transitions to the 
intermediate lifetime regime. Ideally, the results of future searches for stop NLSP will be cast 
in the mass-lifetime plane, yielding a complete characterization of this simple, minimal, and 
well- motivated new physics scenario. 
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A Matrix element for stop decay 

The gravitino field tp^ couples to the supercurrent S'^ as (see, e.g., |l6l HTJ SH] 



J2 iD,4>)^ ^.7 V X + h.c. - [7^ 7l 7^A'^i^;, (A.l) 

^ {'f>,x) 



where (0, x) the scalar and the spinor components of a chiral multiplet and Mp = 
I/VSttG ~ 2.435 X 10^*^ GeV is the reduced Planck mass. The interactions relevant to us 
involve the SU{2) singlet 

= , x = Pnt (A.2) 

and the SU{2) doublet 

Using ci = cos 6i, = sin 6^ to characterize the mixing of the stops, the vertices relevant to 
the lighter stop eigenstate (which would be completely left-handed for = 0) are 

^iwbG = W^t* i^,!"!' CiPLb + h.c. (A.5) 

where we used = + igW^T"- -|- . . . , where g = e/ sin 6'vk- 
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At energies ^ itIq, the gravitino ip^j, can be replaced by the goldstino ip via 

Vv^-,/|^ (A.6) 



Using the gravitino mass expression ttIq = F/y/SMp where \/F is the SUSY breaking scale, 
we get 

Ci,a = J dj* d.i^Yr {ciPL + SiP^) t + h.c. (A.7) 

= ^ 9 W^t* d.i^rr CiPLh + h.c. (A.8) 
We can use {7'^,7'^} = 2g^'^ and = (for an on-shell goldstino) to make the replacement 

di^i^YY = 2 3"^- d^i^-ff'Y 2 d^i^ (A.9) 
which gives (as in eq. (51) of |15] ) 

J^itG = J {ciPL + SiPn) t - d,tt {ciPn + s^P^) ^>^^p] (A.IO) 

and (as in eq. (56) of |18] ) 

^twbG = ^9Ci {W; i* d^^P^b + W- thPa d^^lj) (A.ll) 
Using (lA.7p . (lA.lip . and the standard model interaction 

Ctwb = [W^ tYPLb + W- bYPLt) (A. 12) 
we obtain the matrix element (12. ip quoted in the text. 

B Details relevant to exclusion limits 

B.l ti and tf^ production cross sections 

NNLOapprox+NNLL calculations give the ti production cross section at the Tevatron as ctjj = 
6.30 ± O.lQlolg pb [50] or T.OSloJ^^aP pb [51]. As a compromise between the two sources 
we will assume = 6.7 pb ± 6%. For 7 TeV LHC, the result is a^j = 149 ± 7 ± 8 pb [50] or 
1631^1^ pb [51]. We will use = 156 pb ± 7%. 

NLO+NLL predictions for tt* production cross sections are available in [52]. We model 
their results for the Tevatron by using the leading-order expressions (see, e.g., [53]) with 
MSTW 2008 NLO PDFs multiplied by a i^-factor whose value and uncertainty vary linearly 
between 1.58 ± 13% and 1.48 ± 9% as the stop mass varies between 100 and 200 GeV. For the 
7 TeV LHC we model the results similarly with a K-iactoi of 1.73 ± 13%. The cross sections 
are plotted in figure [H 
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B.2 Systematic uncertainties 

This section describes the systematic uncertainties we assume for the various cut-and-count 
experiments that we analyze. 

For backgrounds other than tt we always use the uncertainties computed by the experimen- 
tal groups. For ti, we also use the result obtained by the experimentalists whenever available. 
In cases where it is not given explicitly we compute it as follows. We include a 6% (Tevatron) 
or 7% (LHC) uncertainty in the theoretical cross section and a 6% (Tevatron) or 3.5% (LHC) 
luminosity uncertainty. For the CDF analyses of [15] (or [16]) we also include 5.0% (or 4.8%) 
acceptance uncertainty (due to the jet energy scale etc.), and for the 6-tagged samples also 
3.5% (4.1%) uncertainty from b tagging efficiency. For the ATLAS analysis [23] we also include 
9% uncertainty from jet energy scale, 7% uncertainty from ISR/FSR, and 11% uncertainty 
from the b tagging efficiency. 

For the ATLAS analysis [20] we take the ti and the other background uncertainties from [20] 
and combine the ti uncertainties in the ee, fifi and e/i channels as correlated and the uncer- 
tainties of the other backgrounds as partly correlated (by taking the average of the correlated 
and uncorrelated result due to lack of more precise information). We consider the uncertainty 
of the other backgrounds to be uncorrelated with that of ti. 

The systematic uncertainty of the stop signal (which is not indicated explicitly in table [3] 
but taken into account even though its effect on the results is very small) is dominated by the 
theoretical uncertainty on the cross section (appendix IB.1|) . We also include the luminosity 
uncertainty quoted above and (in analogy with ti) a 5% (Tevatron) or 10% (LHC) uncertainty 
on the acceptance and 4% (Tevatron) or 10% (LHC) uncertainty from the b tagging efficiency 
where relevant. 

C Simulation details 

For both producing stop pairs and showering the stop decay products, we run Pythia 8.145 [SU 
[55] with Tune 2C [56] and CTEQ 6L1 PDFs [57] for Tevatron analyses and Tune 4C [56] with 
the CTEQ6.5 PDFs [58] (obtained via LHAPDF [59]) for LHC analyses. The intermediate 
step of decaying the stop NLSPs is done according to (12.11) with our own softwarejf] For the 
gravity-mediated stops that we use in the process of deriving limits from the CDF stop mass 
reconstruction procedure of [30] in section 14.21 the decay is performed by Pythia, after the 
desired spectra are created using SuSpect [60] . 

The output of Pythia is further processed with our own software. Lepton identifica- 
tion is simulated by checking whether the leptons pass the various requirements of geometric 
acceptance, px, and calorimeter depositions along their trajectories, as defined in each ex- 
periment. Lepton isolation requirements are also simulated as appropriate in each case. The 
missing energy vector is computed based on the energy deposits in the calorimeters, as 
well as the full energies of muons with |?7| < 2 (CDF or DO) or |?7| < 2.7 (ATLAS). Jets 

^The authors will be willing to provide copies of the software, or parton-level event files produced with it, 
to any interested parties. 
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Analysis 


tt scale factors 


CDF [15], pre-tag sample 
DO [28], up to selection 1 
CDF [30], fe-tagged sample 
ATLAS [2U], pre-tag sample 
ATLAS [H] 


0.88,0.66,0.73 

0.54 
0.90,0.89,0.88 
0.63,1.03,0.74 

1.14 



Table 5: The ratios of the ti event yields predicted by our simulations to those quoted in the 
experimental references, for the various analyses highlighted in section HI Scale factors are 
quoted for ee, fifi and e/i channels separately, where available. 





m-± 




tagless 


tagged 


132.5 


105.8 


47.6 


1.2 


1.3 


135 


125.8 


55.1 


1.2 


2.1 


155.8 


105.8 


64.9 


1.1 


1.1 


160 


125.8 


58.8 


1.1 


1.1 


168.3 


125.8 


71.7 


1.1 


1.2 


179.2 


105.8 


61.3 


0.9 


1.0 



Table 6: Comparison between the results of our and CDF simulation (ratios of our and their 
numbers of events) for gravity- mediated stops studied in [301 EB 132] . 

are clustered using the FastJet [61] implementation of JetClu [62] with cone size R = 0.4, 
overlap threshold 0.75, seed threshold 1 GeV and with ratcheting (for CDF analyses) or the 
anti-^T jet algorithm [63] with cone size R = 0.4 (for ATLAS analyses). We then apply the 
selection requirements, including cuts on px and t] of the various objects, I^t, Ht and other 
analysis-specific quantities. 

For the "tagged" samples, we require at least one h tag, simulated as follows. We determine 
which of the jets that satisfy the analysis conditions are fe-jet candidates by requiring them to 
be within an i? = 0.4 cone of any of the two parton-level h quarks, within ±60% of px of that h 
quark, and have |?7| < 1.5 (CDF) or |?7| < 2.5 (ATLAS). At most two jets are allowed to become 
such candidates. Then we assume fixed tagging efficiency for each of them (we used 0.45 for 
CDF analyses, 0.7 for the ATLAS analysis [20], and 0.5 for the ATLAS analyses [231133]). 

We do not simulate certain factors such as lepton reconstruction efficiencies etc. However, 
we effectively take them into account to some extent by comparing our results for ti with 
those obtained by the experimentalists in their simulations, and multiplying our results for 
the stops by the same scale factors, as described by eq. (14. ip . The scale factors for the analyses 
of table [3] are shown in table [51 Since even the scale factors themselves are not too far from 
1, we estimate that after applying the scale factors our simulated event yields for the stop are 
correct to within ±10%. This is confirmed by table [61 where we compare our results for stops 
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in gravity- mediated scenarios with those obtained by CDF PU| EH [32] (after normahzing our 
predictions to the stop cross sections used by CDF)0 

D CDF stop mass reconstruction algorithm 

We have implemented the dilepton stop mass reconstruction algorithm of [301 Ell 132] as follows: 



The two highest-E'r jets are assigned to the two leptons in the way that results in a 
smaller value of rnj_^f^^ + mj^f^^. 

The possible direction of each Xi + pa-i^ is scanned over the full range. At each 
(01 ! 02) point, the MET carried by each Xi + pa-i^' is then fixed by pr conservation. 

Each Xi + ly pair is modeled as a "pseudoparticle" with mass tuq = 75 GeV, and the 
chargino mass m^± is assumed to be either 105.8 or 125.8 GeV (the values for which 
the results are available in the CDF study). Then the on-shell conditions for Xi ~^ ^ 
and x^ solve the system. This is subject to a 4-fold ambiguity since for each of the two 
possible pairings of (x? + i^) with £ there are two solutions to a quadratic equation. All 
the solutions are kept. 

Find the stop mass m^°^ that minimizes for a given (0i, 02) point, where 

2 / \ 2 



X 



E 

fc=l,2 



^1 



'xt 



(D.l) 



with Fo = 10 GeV, F-i = 2 GeV, Fj = 1.5 GeV. The two pzfi (pseudoparticle z 
momenta) and the two mo are allowed to vary in the fit in order to minimize x^. CDF 
vary also the magnitudes (but not the directions) of the 4-momenta of the leptons and 
jets to account for measurement errors and include the corresponding terms in x^. 



Average the resulting stop mass over all the choices of (0i, 02), weighed by e 



2) 
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